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We have formulated a relaxation mechanism for ferrites and ferromagnetic metals whereby the 
coupling between the magnetic motion and lattice is based purely on continuum arguments con- 
cerning magnetostriction. This theoretical approach contrasts with previous mechanisms based 
on microscopic formulations of spin-phonon interactions employing a discrete lattice. Our model 
explains for the first time the scaling of the intrinsic FMR linewidth with frequency, and temper- 
ature dependence and the anisotropic nature of magnetic relaxation in ordered magnetic materials, 
where M is the magnetization. Without introducing adjustable parameters our model is in reason- 
able quantitative agreement with experimental measurements of the intrinsic magnetic resonance 
linewidths of important class of ordered magnetic materials, insulator or metals. 



PACS numbers: 76.50.-|-g 



INTRODUCTION 



Since the discovery of magnetic resonance, the physics 
community has been fascinated with possible mecha- 
nisms to explain the absorption linewidth or the relax- 
ation time in magnetic materials. It was and still is a very 
challenging problem. Magnetic relaxation is so important 
to understand because it affects a number of technologies, 
including computer, microwave, electronics, nanotechnol- 
ogy, medical, etc.. Ultimately, the physical limitation of 
any technology which incorporates magnetic materials of 
any size, shape and combinations thereof comes down to 
precise knowledge of the relaxation time of the magnetic 
material being utilized. The background of various calcu- 
lations or formulations of magnetic relaxation for the past 
sixty years or so can be summarized briefly as follows: 
(i) The relaxation times in paramagnetic materials [ij 
is characterized by two parameters, Ti and T2, wherein 
describes the magnetic resonance linewidth and Ti 
describes the time taken for the external magnetic field 
Zeemann energy density —Hext • M to relax into thermal 
equilibrium. These times have been modeled in terms 
of various coupling schemes, i.e. spin-spin and/or spin- 
lattice interactions j3l- Since the coupling between spins 
is relatively weak, as it should be in a paramagnetic mate- 
rial, the coupling to the lattice involves discrete spin sites 
rather than a collective cluster of spins. As such, param- 
agnetic coupling is necessarily microscopic in nature. For 
example, a microscopic coupling scheme was formulated 
Q whereby a spin Hamiltonian was modulated by the 
lattice motion. Variants to this approach have been very 
successful in explaining relaxation in paramagnetic ma- 
terials, (ii) The magnetic relaxation of ferrimagnetic or 
ferromagnetic resonance (FMR) linewidth is character- 
ized by the Gilbert parameter [J| a, or equivalently by 
Landau-Lifshitz parameter [3| Al. A distinguishing fea- 



ture of the collective coherent magnetic moments in FMR 
is that the magnitude of the magnetization, M = |M| re- 
mains fixed which requires a magnetic resonance equation 
of the simple form 



— = 7M X Utot 

at 



7M X (H + H'), 



(1) 



wherein the gryomagnetic ratio 7 = ge/2mc. The total 
magnetic intensity Utot has a thermodynamic part H 
determined by the energy per unit volume u, 



du = Tds + H • dM, 



(2) 



and a dissipative part H' determined by the Gilbert lin- 
ear operator a. 



H' 



1 



7M 



(3) 



Eqs.((T|) and ([3]) imply that all components of the magne- 
tization must relax simultaneously in a way which con- 
serves the magnitude of the magnetization. Much of the 
successful microscopic approaches or formulations uti- 
lized in paramagnetic materials were transferred over to 
models loj which attempted to explain Eqs.ll]) and ([3]). 
In some sense this presented a contradiction or paradox 
which was conveniently ignored. As it is well known that 
collective excitations in a ferri or ferromagnetic state can 
be adequately described in classical continuum termi- 
nologies, although microscopic descriptions remain per- 
haps more accurate [vj . To our knowledge very few or any 
microscopic models have been successful in explaining the 
origin of Eq.Q. For example, much attention was given 
in the seventies to explain the FMR linewidth in YIG 
(Y3Fe50i2), since its linewidth was the narrowest ever 
measured in a ferrimagnetic material Clearly, there 
was less to explain, and perhaps spin-lattice interactions 
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could be treated at discrete spin sites as in paramagnetic 
materials. These calculations Q contained many approx- 
imations and predicted an FMR linewidth about 1/10 to 
1/100 of the measured linewidth. We believe that this 
is the best agreement between theory and experiment on 
relaxation in an ordered magnetic material. The pur- 
pose of this work is to improve upon the predictability 
of a theoretical model not only on a given material but 
in general for any ordered magnetic materials without 
restoring to any approximations and assumptions. 

We have adopted a conventional continuum magneto- 
mechanical description of the magnetic and elastic states 
of the ferri or ferromagnetic crystal 0, [l2|. The ad- 
vantage of this description is that the microscopic spin- 
lattice coupling need not be formulated, since it has al- 
ready been included in the continuum model which has 
been proved to be experimentally correct. We introduce a 
thermodynamic argument stating that the heat exchange 
between the magnetic and elastic systems must be the 
same. As such, Eq.(l3]) may be directly related to the elas- 
tic sound wave relaxation time and the coupling strength 
between the magnetic and elastic systems. Specifically, 
we will show that a is proportional to the square of the 
magnetostriction constant, i.e. and inversely propor- 
tional to jMt wherein r the elastic relaxation time. In 
addition, the model predicts that a cannot be presumed 
to be a scalar as it has been done in the past; i.e. a 
is predicted to be anisotropic a second rank tensor in a 
single crystal material. 

It is clear that one needs an interaction between 
phonons and electron spins to account for Gilbert damp- 
ing parameter a. Suhl [lol | and more recently Hickey 
and Moodera [ijj have considered such coupling schemes. 
The Gilbert damping parameter can be thought of as a 
transport coefficient in much the same way as conduc- 
tivity and/or viscosity are transport coefhcients. Such 
transport coefficients describe heating processes by which 
otherwise long lived modes are damped. One can in 
fact relate the Gilbert damping parameter to conduc- 
tivity and/or viscosity. For metallic ferromagnetic mate- 
rials, conductivity as well as electron viscosity produces 
considerable amount of magnetic damping via eddy cur- 
rent heating. For magnetic insulators it is the viscos- 
ity which determines the magnetic damping. As it is 
well known, conductivity and viscosity can be non zero 
even in zero frequency limit. Hence, the implied Gilbert 
damping parameter is also non zero at zero frequency. In 
Suhl and Hickey and Moodera's papers they find, in the 
limit of zero frequency and zero wave number, that the 
real part of a is zero. This limiting case suggests that 
they have not included the zero frequency transport co- 
efficients consistently in their theory. In our derivation 
the expected result at zero frequency occur naturally in 
our formalism. In general, we believe the very nature of 
discreteness (as in paramagnetic materials) gives rise to 
relatively long magnetic relaxation times. However, the 



magnetic relaxation time of a coherent collection of spins 
(as in FMR) implies shorter relaxation times, since it in- 
volves collective acoustic waves in the interaction scheme. 
Our present theoretical treatment takes this into account 
via the continuum magneto-mechanics. 



THEORETICAL MODEL 

From Eq.([3]), it is evident that the heating rate per 
unit volume due to the dissipative magnetic intensity H' 
obeys 



Q 



7M 



M 



dM ^ dM 
"dt IT 
M 



Q — — NiUijNj wherein N = — , 
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and (5 is a second rank tensor 



M 



'■zy 



OLzz 



(4) 



(5) 



The crystal displacement u yields in elasticity theory [ic 
the strain tensor 



(6) 



In virtue of the magneto-elastic effect [l4|, a chang- 
ing magnetization dM-ldt will produce a changing strain 
de/dt. In detail, in terms of third rank magneto-elastic 
tensor Aijki one finds 



= KjkiNkNi, 
Bij = 2XijkiNkNi. 



(7) 



Finally, the fourth rank crystal viscosity tensor, r]ijki 
dtermines the heating rate per unit volume due to the 
time dependent strain 



^ij Vijkl ^kl ■ 



(8) 



Employing Eqs.© and ^ and comparing the result to 
Eq.(|3]) yields the central result of our model. 

For any crystal symmetry the Gilbert damping tensor 
due to magnetostriction coupling is rigorously given by 



47 
M 



(9) 



The following properties of the Gilbert damping tensor 
Eq.® are worthy of note: (i) The Gilbert damping ten- 
sor a is inversely proportional to the magnetization mag- 
nitude M . (ii) The Gilbert damping tensor a is propor- 
tional to the squares of the magnetostriction tensor el- 
ements, (iii) The tensor nature of a dictates that the 
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magnetic relaxation is anisotropic. To a sufBcient degree 
of accuracy, one may employ an average of the form 



frequency / as 



1 



tr{a} 



ax 



(10) 



defining a scalar function a. (iv) The crystal viscosity 
tensor rjnmri may be employed to describe the acoustic 
wave damping [15] . For a mode label a, e.g. a longitudi- 
nal (a = L) or a transverse (a = T) mode, the acoustic 
absorption coefficient at frequency ld is given by [l^ 



-1 _ ^ rja 



(11) 



wherein Va is the acoustic mode velocity and p is the 
mass density. Finally, for a cubic crystal, there are only 
two independent magneto-elastic coefficients which may 
be defined 

3 3 

^xxxx — 2-^100 and A^yxy = 2^^^ ^^'^'^ 

wherein the Cuachy three index magneto-elastic coeffi- 
cients are Xijk ■ 



COMPARISON WITH EXPERIMENT 

The Gilbert damping factor a may be deduced from 
the measurement of the intrinsic FMR linewidth. How- 
ever, the measurement of the intrinsic linewidth is, in- 
deed, very difhcult. The reason for this conclusion is 
that there are too many extrinsic effects that influence 
the measurement. For example, in ferromagnetic metals 
like Ni, Co and Fe the intrinsic linewidth contribution to 
the total linewidth measurement 3, 17 1 may be between 
10% and 30%. The rest of the hnewidth fl^ may be due 
to exchange-conductivity effects. 

However, there may be other contributions, such as 
magnetostatic excitations, surface roughness, volume de- 
fects |19|, crystal quality, interfaces [20|, size, etc.. Sim- 
ilar conclusions apply to ferrites except there are no 
exchange-conductivity effects [Si- Thus, the reader 
should be mindful that when we quote or cite an intrinsic 
value of the linewidth it represents a maximum value for 
there can be some hidden extrinsic contributions in an 
experiment. However, we have relied on data well estab- 
lished over the years. The criteria that we have adopted 
in choosing an ensemble of intrinsic linewidth measure- 
ments are the ones exhibiting the narrowest linewidth 
ever measured in single crystal materials. In addition, 
we required full knowledge of their elastic, magnetic and 
electrical properties [H, [l3, El, lllj . The objective is not 
to introduce any adjustable parameters. 

The experimental value of Gilbert damping parameter 



(13) 



The factor ^3/2 assumes Lorentzian line shape of the res- 
onance absorption curve. The theoretical Gilbert damp- 
ing parameter ath value is expressed in terms of known 
[l7| parameters so that there are no adjustable parame- 
ters in our comparison to experiments, as shown in TA- 
BLE IT] The theoretical prediction for the Gilbert damp- 
ing paramter is that 



ath 



36p7 
Mr 



>2 



^111 



9 

9t 



(14) 



ttcxp may be deduced from the FMR linewidth AiJ at 



wherein p is the mass density, qr ~ ^Tj^ is the trans- 
verse acoustic propagation constant, is the longitudi- 
nal acoustic propagation constant, vr is the transverse 
sound velocity, A is the exchange stiffness constant, Aioo 
and Am are magnetostriction constants for a cubic crys- 
tal magnetic material. The transverse acoustic propaga- 
tion constant, was approximated on the basis that the re- 
laxation process conserved energy and wave vector. Since 
the acoustic frequency is fixed in the process the longi- 
tudinal propagation constant may be also calculated to 
be qL = qri'^T/^L) for magnetic materials, wherein vl 
is the longitudinal sound wave velocity. 

In FIG in we plot the experimental and theoretical val- 
ues Gilbert damping constants as given by Eas. ([T^ and 
p4| . We note that the agreement between theory and 
experiment is remarkable in view of the fact that any 
of the cited parameters could differ from the ones listed 
in TABLE U by as much as 20-30%. For example, the 
linewidth reported in TABLE |T] may not be on the same 
sample where the elastic or magnetic parameters were 
cited. In a few cases we needed to extrapolate the value 
of A, since there was no published value. In FIGlU we 
did not present data on the ferromagnetic metals for lack 
of confidence on the linewidth data. For example, mag- 
netostatic mode excitations have a deleterious effect on 
the dependence of the FMR linewidth on size. Most, if 
not all, previous FMR linewidth measurements have been 
performed on slabs, wiskers, etc.. whcih can indeed sup- 
port magnetostatic mode excitations. Additional compli- 
cations arise as a result of exchange-conductivity excita- 
tions in the linewidth data. Nevertheless, the agreement 
between theory and experiment is quite satisfctory. 



CONCLUSION 

Qualitative and quantitatively our model is in agree- 
ment with experimental observations of the intrisic FMR 
linewidth reported over the years. Speciafically, experi- 
mentally the most important characteristics of the intrin- 
sic FMR linewidth, A_ff , measured on ordered magnetic 
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TABLE I: Calculated and measured Gilbert damping {a) parameters 





QT 


Aioo 


Am 


M A 


AH 


/ 


r 


Ctth 




Materials 


(lO^^cm^^) 


(io-«) 


(io-«) 


(G/47r) (lO-'^erg/cm) 


(Oe) 


(GHz) 


(10-i^sec) 


(10-5) 


(10-5) 



Y3Fe50i2 " 


3.8 


1.25 


2.8 


139 


0.40 


0.33 


9.53 


4.4 


5.56 


9.0 


Y3Fe4GaOi2 " 


1.46 


- 1 


- 1 


36 


0.28 


3.0 


9.53 


4.4 


51 


76 


Lio.5Fe2.5O4 ' 


8.6 


- 8 


+ 


310 


0.40 


2.0 


9.50 


1.5 


26 


50 


NiFe2 04 


7.49 


- 63 


- 26 


270 


0.40 


35 


24.0 


710 


26 


350 


MgFe2 04 ' 


9.30 


- 10 


- 1 


90 


0.1 


2.3 


4.9 


1.5 


120 


120 


MnFe2 04 * 


6.6 


- 30 


- 5 


220 


0.4 


238 


9.2 


1.5 


930 


1040 


BaFei2 0i9 " 


9.6 




15 


350 


0.4 


6 


55 


1.5 


18 


26 


Ni'' 


6.3 


- 46 


25 


484 


0.75 


102 


9.53 


1.8 


770 


2600 


Fe"^ 


8.75 


20 


- 20 


1690 


1.9 


9 


9.53 


1.8 


30 


220 


00"^ 


5.1 




80 


1400 


2.78 


15 


9.53 


1.8 


530 


370 



"Garnets ''Spinels ^Hexagonal Ferrite ''Ferromagnetic Materials 
(Note: Longitudinal acoustic wave constant is ql ~ {vt /vL)qT) 



10' 



10"' 



10" 



10 



10"' 



of samples seemed to have incited the controversy. Im- 
proved or more accurate angular Hnewidth data 18, l3| 
supports the notion of an anisotropic Hnewidth in or- 
dered magnetic materials in agreement with our model. 
In summary, we believe that the comparison between 
theory and experiment is very encouraging in terms of 
continuing this continuum approach to explain intrinsic 
linewidths in ordered magnetic materials. 
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FIG. 1: Shown are the experimental and theoretical values of 
the Gilbert damping constants as given by Eas. (|13p and (|14|l . 



materials (metal or insulator) for the past fifty years are 
that AT? scales with frequency and [3, IHj HI] ■ In- 
deed, these are the predictions of our theory. In addition, 
AT? scales with the magnetostriction constant squared, 
see FIGlH FIG H] was plotted in a logarithmic scale only 
to be able to include all of the data in TABLE [H Another 
prediction of our theoretical work is that the Gilbert 
damping parameter a is not simply a scalar parame- 
ter but a tensor quantity. This implies that the FMR 
hnewidth is intrinsically anisotropic in single crystals of 
ferri- ferromagnetic materials. There was much contro- 
versy in the seventies about whether or not the intrin- 
sic linewidth should be anisotropic or not. Poor quality 
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